SPhT-T05/241 



Formal matrix integrals and combinatorics of maps 



B. EynarJ^ Service de Physique Theorique de Saclay, 



O 

o 

Csl ; F-91191 Gif-sur-Yvette Cedex, France. 

O 
D 

Q 



Abstract: 

This article is a short review on the relationship between convergent matrix inte- 
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^ : 1 Introduction 

^ I This article is a short review on the relationship between convergent matrix inte- 

^ ', grals, formal matrix integrals, and combinatorics of maps. We briefly summa- 

rize results developed over the last 30 years, as well as more recent discoveries. 

We recall that formal matrix integrals are identical to combinatorial generating 
^ ' functions for maps, and that formal matrix integrals are in general very different from 

d I convergent matrix integrals. Both may coincide perturbatively (i.e. up to terms smaller 

than any negative power of A^) , only for some potentials which correspond to negative 
weights for the maps, and therefore not very interesting from the combinatorics point 
of view. 

We also recall that both convergent and formal matrix integrals are solutions of the 
same set of loop equations, and that loop equations do not have a unique solution in 
general. 

Finally, we give a list of the classical matrix models which have played an important 
role in physics in the past decades. Some of them are now well understood, some are 
still difficult challenges. 

Matrix integrals were first introduced by physicists [SS], mostly in two ways: 
- in nuclear physics, solid state physics, quantum chaos, convergent matrix integrals are 
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studied for the eigenvalues statistical properties [HI [331 [3 ISS]. Statistical properties 
of the spectrum of large random matrices show some amazing universal behaviours, 
and it is believed that they correspond to some kind of "central limit theorem" for 
non independent random variables. This domain is very active and rich, and many 
important recent progresses have been achieved by the mathematicians community. 
Universality was proved in many cases, in particular using the Riemann-Hilbert ap- 
proach of Bleher-Its [6] and Deift Venakides Zhou Mac Laughlin [19], and also by large 
deviation methods [SI [35] . 

- in Quantum Chromodynamics, quantum gravity, string theory, conformal field theory, 
formal matrix integrals are studied for their combinatorial property of being generat- 
ing functions of maps [20]. This fact was first discovered by t'Hooft in 1974 [19], 
then further developed mostly by BIPZ ^12j as well as Ambjorn, David, Kazakov 
[201 HSl E21 [371 [3B]. For a iong time, physicist's papers have been ambiguous about the 
relationship between formal and convergent matrix integrals, and many people have 
been confused by those ill-defined matrix integrals. However, if one uses the word "for- 
mal matrix integral", many physicist's results of the 80's till now are perfectly rigorous, 
especially those using loop equations. Only results regarding convergency properties 
were non rigorous, but as far as combinatorics is concerned, convergency is not an 
issue. 

The ambiguity in physicist's ill-defined matrix integrals started to become obvious 
when E. Kanzieper and V. Freilikher [12], and later Brezin and Deo in 1998 [TT] tried 
to compare the topological expansion of a formal matrix integral derived from loop 
equations, and the asymptotics of the convergent integral found with the help of or- 
thogonal polynomials. The two results did not match. The orthogonal polynomial's 
method showed clearly that the convergent matrix integrals had no large N power 
series expansion (it contained some (—1)^). The origin of this puzzle has now been 
understood [9], and it comes from the fact that formal matrix integrals and convergent 
matrix integrals are different objects in general. 

This short review is only about combinatoric properties of formal matrix integrals. 
Matrix models is a very vast topic, and many important applications, methods and 
points of view are not discussed here. In particular, critical limits (which include 
asymptotics of combinatoric properties of maps), the link with integrable systems, 
with conformal field theory, with algebraic geometry, with orthogonal polynomials, 
group theory, number theory, probabilities and many other aspects, are far beyond the 
scope of such a short review. 
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2 Formal matrix integrals 

In this section we introduce the notion of formal matrix integrals of the form: 

^wai(i) = / dM, dM,... dM, fe^i % ^ -NViM„...,M.)) (2-1) 

J formal 

The idea is to formally expand the exponential in powers of t~^, and compute 

the Gaussian integral for each term. The result is a formal series in powers of t. So, 
let us define it precisely. 



Definition 2.1 Q is an invariant non- commutative monomial of Mi, . . . , Mp, ifQ — 1 
or if Q is of the form: 

- ^ 1 

^ = Un^' (^r) (2-2) 

r=l 

where each Wr is an arbitrary word written with the alphabet Mi, . . . ,Mp. Q is the 
equivalence class of Q under permutations of the Wr 's, and cyclic permutations of the 
letters of each Wr- 

The degree of Q is the sum of lengths of all Wr 's. 

Invariant non- commutative polynomials of Mi,...,Mp are complex finite linear 
combinations of monomials: 

V = Y,tQQ , tQ^C (2-3) 
Q 

The degree of a polynomial is the maximum degree of its monomials. 

They are called invariant, because they are left unchanged if one conjugates all 
matrices Mj UMjU~^ with the same invertible matrix U . 

Invariant polynomials form an algebra over C. 

Let V{Mi, . . . , Mp) be an arbitrary invariant polynomial of degree d in Mi, . . . , Mp, 
which contains only monomials of degree at least 3. 

Proposition 2.1 LetC be apxp symmetric positive definite matrix, then the following 
Gaussian integral 

where dMi is the usual Lebesgue (U (iV) invariant) measure on the space of hermitian 
matrices H^, is absolutely convergent and has the following properties: 
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Ak{t) is a polynomial in t, of the form: 



Akit) 



(2-5) 



k/2<j<kd/2-k 




proof: 



y4o = 1 is trivial. Let d = deg V. Since V is made of monomials of degree at least 3 
and at most d, then is a sum of invariant monomials whose degree / is between 3k 
and dk. According to Wick's theorem, the Gaussian integral of a monomial of degree I 
is zero if I is odd, and it is proportional to t'^^ if / is even. Since 3k < I < dk we have: 



Thus Ak{t) is a finite linear combination of positive integer powers of t, i.e. it is a 
polynomial in t, of the form of eq J2-5l 

The matrix size A^'s dependence comes in several ways. First there is the factor 
A^^'^. The matrix size also appears in the matrix products (each matrix product is a 
sum over an index which runs from 1 to A^), in the traces (it means the first and last 
index of a matrix product have to be identified, thus there is a Kroenecker's 6ij of two 
indices). And after Gaussian integration over all matrix elements, the Wick's theorem 
pairings result in A^~'/^ times some product of Kroenecker's S of pairs of indices (times 
some elements of the matrix which are independent of A^). The matrix indices 
thus appear only in sums and 6's, and the result of the sum over indices is an integer 
power of A^. Thus, each A^it) is a finite sum (sum for each monomial of V'', and the 
Gaussian integral of each monomial is a finite sum of Wick's pairings) of positive or 
negative powers of A^, i.e. a Laurent polynomial in A^. □ 

Definition 2.2 The formal matrix integral ^formai(^) is defined as the formal power 



< k/2 <l/2 ~ k < dk/2 - k 



(2-6) 



series: 





(2-7) 



fe=0 



and each Zi is a Laurent polynomial in N. Notice that Zq = 1. 



By abuse of notation, Zformaiit) is often written: 



but it does not mean that it has anything to do with the corresponding convergent (if 
it converges) integral. In fact, the integral can be absolutely convergent only if deg(l^) 
is even and if the tg corresponding to the highest degree terms of V have a negative 
real part. But as we shall see below, the relevant case for combinatorics, corresponds 
to all tg's positive, and in that case, the formal integral is NEVER a convergent one. 

Definition 2.3 The formal free energy Fformai(i) is defined as the formal log of Z formal- 



We have Fq = 0. Each Fj is a Laurent polynomial in N. 
2.1 Combinatorics of maps 

Recall that an invariant monomial is a product of terms, each term being the trace of 
a word in an alphabet of p letters. Thus, an invariant monomial is given by: 

• the number R of traces, (i? — 1 is called the crossing number of Q), 

• R words written in an alphabet of p letters. 

The R words can be permuted together, and in each word the letters can be cycli- 
cally permuted. We label the invariant monomials by the equivalence classes of those 
permutations. 

Another graphical way of coding invariant monomials is the following: 

Definition 2.4 To each invariant monomial Q we associate biunivoquely a Q — gon 
(generalized polygon) as follows: 

• to each word we associate an oriented polygon (in the usual sense), with as many 
edges as the length of the word, and whose edges carry a "color" between 1 andp, given 
by the corresponding letter in the word. 

• the R words are glued together by their centers on their upper face ( in accordance 
with the orientation) , so that they make a surface with R — 1 crossings. 

• R — 1 which is the number of traces minus one (i.e. one trace corresponds to a 
crossing number zero), is called the crossing number of the Q-gon. 

• The degree deg{Q) of the Q-gon is the total number of edges (sum of lengths of all 



• to Q-gon we associate a symmetry factor sq = #Aut(Q) which is the number of 
symmetries which leave Q unchanged. 

An example is given in figHJ Notice that we allow a Q-gon to be made of polygons 
with possibly one or two sides. We will most often call the Q-gons polygons. The usual 
polygons are Q-gons with no crossing i.e. R = 1. 




(2-9) 



words). 
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Figure 1: The invariant monomial Q = N'^ Tr(Aff Afa) Tr(A'/|) Tr(Af|Af4Aff ) of degree 
11, is represented by 2 triangles and one pentagon glued together by their center. The 
dashed lines mean that the 3 centers should actually be seen as only one point. Its 
symmetry factor is sq = 3 because we can perform 3 rotations on the triangle of color 
(3,3,3). 



Definition 2.5 Let p > 1 and d > 3 be given integers. Let Q^^p be the set of all 
invariant monomials Q (or polygons) of degree 3 < deg{Q) < d. 
Qd,p is clearly a finite set. 

Definition 2.6 Let Si^d,p be the set of oriented discrete surfaces such that #edges- 
#Qgons+#crossings = I, and obtained by gluing together polygons (belonging to Qd,p) 
by their sides (following the orientation) . The edges of the polygons carry colors among 
p possible colors (thus each edge of the surface, is at the border of two polygons, and 
has two colors, one on each side). 

Let Si^d-p be the subset of Si^d,p which contains only connected surfaces. 

Such surfaces are also called "maps". 

Proposition 2.2 Si^d,p is a finite set. 
proof: 

Indeed, since all polygons of Qd,p have at least 3 sides, we have #edges > 
|#polygons, and thus #edges < 21 and ^polygons < 4Z/3, and thus the number 
of discrete surfaces in Si, is finite. We can give a very large bound on ^Si^d,p- We have: 

#Si,d,p < {Ml/3f (2d//3)! {4dl/3)\ < {Adl/Sf {2dl)\ (2-10) 

□ 
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Theorem 2.1 (t'Hooft 1974 and followers): If the potential V is an invariant polyno- 
mial given by 

v= J2 —Q (2-11) 

then: 

j: ^A^Arx(^^nc^'^n((^-%)'"'^^^ (2-13) 

where Zi and Fi were defined in def \2.^ def \2.3[ and where: 

• UqIS) is the number of Q-gons in S, 

• Eij{S) is the number of edges of S which glue two polygons whose sides have colors 
i and j , 

• x{S) is the Euler characteristic of S. 

• Aut(5') is the group of automorphisms of S, i.e. the group of symmetries which 
leave S unchanged. #Aut(S') is called the symmetry factor. 

In other words, Zi is the formal generating function which enumerates discrete sur- 
faces of Si^d,p, according to their Euler characteristics, their number of polygons, num- 
ber of edges according to their color, number of crossings... Fi is the formal generating 
function for the same surfaces with the additional condition that they are connected. 

An example is given in fig. [2l 
proof: 

It is a mere application of Feynman's graphical representation of Wick's theorem^]. 

The only non-trivial part is the power of N, because N enters not only in the 
weights but also in the variables of integrations. It was the great discovery of t'Hooft 
to recognize that the power of is precisely the Euler characteristics. Let us explain 
it again. 

First, V'' is decomposed into a combination of monomials: 

v = y2-Q , v' = y2n,GG (2-14) 

where G is a product of Q's, and Tk,G = T.Q,uQ2U...uQk=G tQitQ2 ■■■tQk n^^- ^ ^^^^ 
a collection of polygons, some of them glued by their centers. So far the polygons of 



^Although Feynman's graphs are sometimes regarded as non-rigorous, let us emphasize that it 
is only when Feynman's graphs and Wick's theorem are used for functional integrals that they are 
non-rigorous. Here we have finite dimensional Gaussian integrals, and Wick's theorem and Feynman's 
graphs are perfectly rigorous. 
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Figure 2: UV = TiMf + |^ Tr M| + TiMff, the above surface contributes 
to the term )f:4 ^4 ^4,4 (C~"^)f (C~^)^ 2 • Indeed Mi is represented in blue, M2 in 
red, so that Tr Mf corresponds to blue squares, Tr M| corresponds to red squares, 
and {Tt Mfy corresponds to pairs of blue squares glued at their center. Its Euler 
characteristic is x = (it is a torus), and this surface has no automorphism (other 
than identity), i.e. #Aut=l. It corresponds to / = 7. 



G are not yet glued by their sides to form a surface. Remember that each Q carries a 
factor A^~^ if Q is made of R traces. 

Then, for each G, the Gaussian integral is computed by Wick's theorem and gives 
a sum of Wick's pairings, i.e. it is the sum over all possible ways of pairing two M's, 
i.e. it is the sum over all possible ways of gluing together all polygons by their sides, 
i.e. corresponds to the sum over all surfaces S. The pairing < MiabMjcd > of the (a, b) 
matrix element of Mj and the (c, d) matrix element of Mj gives a factor: 

< MiabMj,d >= ^ {C-')^, 5adhc (2-15) 

The double line notation for pairings (see [20]) allows to see clearly that the sum 
over all indices is JV#vertices(5)^ The total power of is thus: 

2k — Rq + ^vertices — #edges — ^polygons 
Q 

Now, notice that 

^polygons = ^^Rq = k + ^(-Rq - 1) = A; + #crossings (2-17) 
Q Q 



(2-16) 
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Thus the total power of N is: 

^vertices — #edges + ^polygons — 27^crossings = x (2-18) 

which is the Euler characteristic of S. 

We leave to the reader to see in the hterature how to find the symmetry factor. □ 

Corolleiry 2.1 N~'^Fi is a polynomial in N~'^: 

Fi = Y, Fi,9 (2-19) 

5=0 

proof: 

Indeed the Euler characteristic of connected graphs is of the form x — '^ — '^g where 
g is positive integer and is the genus (number of handles). The fact that it is a 
polynomial comes from the fact that F; is a finite sum. 

It is easy to find some bound on g. We have: 

2g — 2 — — ^vertices + #edges — ^polygons + 2#crossings 
= — ^vertices + 21 — hedges + ^polygons 
(2 - 20) 

and using #edges > 3/2#polygons, we have: 

2^ - 2 < -^vertices + 21- ^#polygons < 21 (2-21) 

□ 

2.2 Topological expansion 

Definition 2.7 We define the genus g free energy as the formal series in t: 

oo 

F(^) {t)^-Y^ Fi,g (2-22) 

1=0 

F^^^ is the generating function for connected oriented discrete surfaces of genus g. 

Remark: the minus sign in front of F^^^ is there for historical reasons, because in 
thermodynamics the free energy is — In Z. 

There is no reason a priori to believe that F^^^t) might have a finite radius of 
convergence in t. However, for many classical matrix models, it is proved that \/g, F^^^ 
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has a finite radius of convergence because it can be expressed in terms of algebraic 
functions. 

There is also no reason a priori to believe that the F^^^^s should be the limits of 
convergent matrix integrals. There are many works which prove that the number of 
terms in F^^^^ grows with g like {Pg) \ for some p. If t and all tg's are positive (this is the 
interesting case for combinatorics and statistical physics because we want all surfaces 
to be counted with a positive weight), then F^^^ is positive and grows like {Pg^-, and 
therefore the sum of the topological series CANNOT be convergent (even under Borel 
resummation). For negative t, it is only an asymptotic series, and at least in some 
cases, it can be made convergent using Borel transforms. 



3 Loop equations 

Loop equations is the name given to Schwinger- Dyson equations in the context of 
random matrices [53j. The reason for that name is explained below. Let us recall that 
Schwinger-Dyson equations are in fact nothing but integration by parts, or the fact 
that the integral of a total derivative vanishes. 

In particular, in the context of convergent matrix integrals we have: 

= V /" dMidM2 ... dMr, I T^T-^T i^^r-^ — I 

((G(Mi, . . . , M,) + G{M„ M,)% ^ -w(m„...,m,))^ 

+ V / dMidM2... dMp — — 

(^(G(Mi, . . . , M,) + G(Mi, . . . , M,)% e-^te- ^ ^^^^ -w(m„...,m,))^ 

(3-1) 

where G is any non commutative polynomial, and k is an integer between 1 and p. 

Therefore, Schwinger-Dyson equations for matrix integrals give relationships be- 
tween expectation values of invariant polynomials. Namely: 

(jll ^^^^^ - ND,{V)^ G^^ = {Ku{G)) (3-2) 

where Dk is the non commutative derivative with respect to Mfc, and Kk{G) is some 
invariant polynomial which can be computed by the following rules: 
• Leibnitz rule 

Kk{A{Mi, . . . , Mfc, . . . , Mp)5(Mi, . . . , Mfc, . . . , Mp)) 
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= {Kk{A{Mi, ...,Mk,..., Mp)B{Mu . . . , m^, . . . , Mp)))^^^^,,^ 

+ {Kk{A{Mu ...,mk,..., Mp)B{Mu . . . , M^, . . . , M^)))^^^^,^ (3-3) 

• split rule 

Kk{A{Mi, . . . , mfc, . . . , Mp)MlB{M,, . . . , m^, . . . , Mp)) 
i-i 

= J] Tr (A(Mi, . . . , m,, . . . , Mp)Ml) Tr (M^-^-^5(Mi, . . . , m,, . . . , Mp))(3-4) 
i=o 

• merge rule 

irfc(A(Mi, . . . , rrtfc, . . . , Mp) Tr {MIB{M,, . . . , m^, . . . , M^))) 
z-i 

= J] Tr (A(Mi, . . . , m,, . . . , Mp)MlB{M,, . . . , m^, . . . , Mp)M^-^-^) (3-5) 

• no Mfc rule 

irfc(A(Mi,...,mfc,...,Mp)) = (3-6) 

Since each rule allows to strictly decrease the partial degree in M^, this set of rules 
allows to compute Kk{G) for any G. 

For any G and any k we get one loop equation of the form eq j3-2[ 

Definition 3.1 The formal expectation value of some invariant polynomial 
P{Mi, . . . , Mk) is the formal power series in t defined by: 

Ak,p[t) - Nry^^P r M M 

Ih.x...xh. dMi dM,... dMp e-^^ ^..^.c,m.m, 
Ak,p{t) is a polynomial in t, of the form: 

AkAt)= ^^^P'^^' (3-8) 

deg(P) /2+fc/2<j<deg(P) /2+kd/2-k 

and we define the following quantity 

2j-dcgP 
k=0 

and the formal series 

oo 

Ap{t)= Apjt^ (3-10) 

j=degP/2 

Again, each A^^p, A^^pj, Apj is a Laurent polynomial in N. 
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The formal expectation value of P is defined as: 



(P(Mi,...,M,)),_, = -M^- (3-11) 



where the division by Zformai is to be taken in the sense of formal series, and it can be 
performed since ^formal (^) = Yl^=o^jt'' with Zq = 1. 

The formal expectation value is often written by abuse of notations 

i , dMr dM,... dM, P e-Tte^i ^TrM.M, -w(m.,...,m,)) 

< -r >formal— TT? ~ ^ — W'^^) 

j,^^^^,dM,dM, . . . dM, e-Tfe-^TrM.M,-W(M„...,M,)) 

Theorem 3.1 The formal expectation values of the formal matrix integral satisfy the 
same loop equations as the convergent matrix integral ones, i.e. they satisfy eq \3-l!\ for 
any k and G. 

proof: 

It is clear that Gaussian integrals, and thus formal integrals satisfy eqs J3-l[ The re- 
mainder of the derivation of loop equations for convergent integrals is purely algebraic, 
and thus it holds for both convergent and formal integrals. □ 

On a combinatoric point of view, loop equations are the generalisation of Tutte's 
equations for counting planar maps [501 El] • This is where the name "loop equations" 



comes from: indeed, similarly to theorem J2.1l formal expectation values of traces are 
generating functions for open discrete surfaces with as many boundaries as the num- 
ber of traces (and again the power of is the Euler characteristic of the surface). 
The boundaries are "loops" , and the combinatorical interpretation of Schwinger-Dyson 
equations is a recursion on the size of the boundaries, i.e. how to build discrete surfaces 
by gluing loops a la Tutte [501 EI] ■ 

Notice that in general, the loop equations don't have a unique solution. One 
can find a unique solution only with additional constraints not contained in the loop 
equations themselves. Thus, the fact that both convergent and formal matrix models 
obey the same loop equations does not mean they have to coincide. Many explicit 
examples where both do not coincide have been found in the literature. It is easy to 
see on a very simple example that Schwinger-Dyson equations can't have a unique 
solution: consider the Airy function j^eh*^~^^ dt where 7 is any path in the complex 
plane, going from cxd to 00 such that the integral is convergent. There are only two 
homologicaly independent choices of 7 (one going from +00 to e^*'^/'^oo and one from 
+00 to e~^*'^/^oo). Schwinger-Dyson equations are: < nt"'~^ + 1"^^ — xt" >= for all 
n. It is clear that loop equations are independent of the path 7, while their solution 
clearly depends on 7. 
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Theorem 3.2 The formal matrix integral is the solution of loop equations with the 
following additional requirements: 

- the expectation value of every monomial invariant is a formal power series in 

- the t ^ limit of the expectation value of any monomial invariant of degree > 1 
vanishes: 

lim< g >= if Q ^ 1 (3-13) 

proof: 

The fact that all expectation values have a formal N""^ expansion follows from the 
construction. The fact that limj_»o < Q >= if Q 7^ follows from the fact that 
we are making a formal Taylor expansion at small t, near the minimum of the quadratic 
part, i.e. near Mj = 0, i = 1, . . . ,p. □ 

However, even those requirements don't necessarily provide a unique solution to 
loop equations. Notice that there exist formal solutions of loop equations which satisfy 
the first point (there is a formal A^~^ expansion), but not the second point (limt^o < 
Q >= 0). Those solutions are related to so-called "multicut" solutions, they also have 
a known combinatoric interpretation, but we don't consider them in this short review 
(see for instance [SH [9] for examples). 

There is a conjecture about the relationship between different solutions of loop 
equations: 

Conjecture 3.1 The convergent matrix integral (we assume V to be such that the 
integral is convergent, i.e. the highest tq 's have negative real part) 

Z_ = [ dM, dM,... dM, ^ M.M, -w(M„...,M,)) ^3_^4^ 

I HP 

is a finite linear combination of convergent formal solutions of loop equations (i.e. a 
formal solution of loop equations \nZ = — Y^'^^^ N'^~'^^ F^^\ such that the N'"^ series 
is convergent.), i.e. 

oo 

^conv = 5ZQZ, , InZ, = -^iV2-29F/^) (3-15) 

i g=0 

Hint: It amounts to exchanging the large N and small t limits. First, notice 
that convergent matrix integrals are usualy defined on if^, but can be defined on any 
"integration path" in the complexified of if^, which is Mjv(C)^, as long as the integral 
is convergent. The homology space of such contours is of finite dimension (because 
there are a finite number of variables p x iV^, and a finite number of possible sectors 
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at oo because the integration measure is the exponential of a polynomial). Thus, the 
set of "generalized convergent matrix integrals" defined on arbitrary paths, is a finite 
dimensional vector space which we note: Gen. The hermitian matrix integral defined 
on is only one point in that vector space. 

Second, notice that every such generalized convergent matrix integral satisfies the 
same set of loop equations, and that loop equations of type eq J3-ll are clearly linear in 
Z. Thus, the set of solutions of loop equations is a vector space which contains the 
vector space of generalized convergent matrix integrals. 

Third, notice that formal integrals are solutions of loop equations, and therefore, 
formal integrals which are also convergent, belong to the vector space of generalized 
convergent matrix integrals. 

Fourth, one can try to compute any generalized convergent matrix integral by small 
t saddle point approximation (at finite A^). In that purpose, we need to identify the 
small t saddle points, i.e. a set of matrices (M^, . . . , Mp) G M^i^Cy such that k 
one has: 

OMT =0 (3-16) 

''''' ^''"^ ^ M„=M„ 

and such that 

Im I 5^ Tr {^Ci,mMiMj - V{M,, . . . , M^) j = Im | ^ Tr {^Q^mK^Mj - V{M„ ...,Mp) 

\ l,m / \ l,m 



(3-17) 



and 
Re 



( {\Ci,^MiM^) - \/(Mi, . . . , Mp) j > Re f 5^ Tr (^G,^M,M„) - V{M^, . . . , M^) 

\ l,im / \ l,m 



(3-18) 

If such a saddle point exist, then it is possible to replace exp hy its Taylor series 
expansion in the integral and exchange the summation and the integration (because 
both the series and the integral are absolutely convergent, this is nothing but WKB 
expansion). This proves that saddle points are formal integrals and at the same time 
they are generalized convergent integrals, thus they are formal convergent integrals. 

The conjecture thus amounts to claim that saddle points exist, and that there exist 
as many saddle points as the dimension of Gen, and that they are independent, so that 
the saddle points form a basis of Gen. 

Notice that a linear combination of convergent formal solutions has no A^~^ expan- 
sion in general, and thus the set of convergent formal integrals is not a vector space. 
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This conjecture is proved for the 1-matrix model with real potential [HI [22], and for 
complex potentials it can be derived from Bertola-Man Yue [3] (indeed, the asymptotics 
of the partition function follow from those of the orthogonal polynomials). It is the 
physical intuition that it should hold for more general cases. It somehow corresponds 
to the small t saddle point method. Each saddle point has a WKB expansion, i.e. is 
a convergent formal solution, and the whole function is a sum over all saddle points. 
The coefficients of the linear combination reflect the homology class of the path (here 
Hff) on which Z is defined. This path has to be decomposed as a linear combination 
of steepest descent paths. The coefficients of that linear combination make the q's of 
eq J3-15l (this is the generalisation of the idea introduced in [9j). 



4 Examples 

4.1 1-matrix model 

where V(M) = f - V{M) is a polynomial in M, such that V{M) contains only 
terms with degree > 3. 

For any a and 7, let us parametrise the potential V as: 



x{z) = a + -y{z + 1/z) 
V'(x(^)) = E-=o^.(^^' + ^"') 



d , (4-2) 



We determine a and 7 by the following conditions: 

vo = , vi = - (4-3) 

7 

i.e. a and /3 are solutions of an algebraic equation and exist for almost any t, and they 
are algebraic functions of t. In general, they have a finite radius of convergence in t. 
We chose the branches of a{t) and 7(t) which vanish at t = 0: 

a(t = 0) = , 7(t = 0) = (4-4) 

Then we define: 



The curve {x{z), y{z)) z E C is called the spectral curve. It is a genus zero hyperellip- 
tical curve = Polynomial(x). It has only two branch points solutions of x'{z) = 
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which correspond to z = ±1, i.e. x = a ± 27. y as a function of x has a cut along the 
segment [a — 27, a + 27]. Notice that we have: 

Res ydx = t = — Res ydx (4-6) 
Res V'{x)ydx = , Res xV'{x)ydx = (4-7) 

2— >00 2— >00 

Then one has 0: 

= ^ f Res V(x)ydx - t Res V(x)— - 1^ -tHn(^ 

2 y2^oo 2^00 Z 2 V ^ 

= E T^^^-+^ - - ^(-l)^(-2,-i - v,,^,) - - In (Cf/t) 

(4-8) 

and: 



ir(i) = __Lin(l^lifAz£Z] (4-9) 



lj^/^7V(l)2/'(-l) 
24 V ^' 

Expressions are known for all F^^^^s and we refer the reader to [231 E]- Those 
expressions are detailed in the next section about the 2-matrix model, and one has 
to keep in mind that the 1-matrix model is a special case of the 2-matrix model with 

Vi(x)=V(x) + 4andV2(2/) = 4- 

As a fore-taste, let us just show an expression for F^^^: 

_2f{2) = Res Res Res $(21)^(21,22)^^(1/^1,^3^ ^ ^ 



±\2 (^^ _ 1\2 



(4 - 10) 



+2 Res Res Res $(21)^(21,22)^(^2,^3)71 r\?7 TT? 

2l^±l 22^±l 23-+±l -i i-)'!' [Zn -V 

\ 21 -j^ ^2 ^ 1 

+2 Res Res Res $(21)^(21,22)^(22,23)71 ^71 

2i^±i 22^±i 23^±i H 23 r 

^21 23' ^22 



where the residue is first evaluated for 23 then 22 then 21, and where: 

"'^ ^ 4^ ^^(T^TkT^^m^ ^^"^^^ 

$(2) = — -r- I ydx (4-12) 



47 zy{z){z - i) ii/ 

4.1.1 Example triangulated maps 

Consider the particular case V{x) = ^x^ 
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Let T — and a be a solution of: 

a-o? = AT (4-13) 
and consider the branch of a{T) which goes to 1 when T = 0. If we parametarize: 

T = -| = ^ 4-14 

C3 6^3 



we have: 



We have: 



which gives: 



cos f — , , 

a = 4-15 

cos (f ) ^ > 

«-|(l-a) , (4-16) 

^^0 = , vi = - , -^37^ (4-17) 

7 



^U)_JLi,( ^-^"V^-' ] (4.19) 



1 , /l-2a\/l-a^ 

— m ; 

24 V « 

The radius of convergence of F^^^ is \T\ < for all g. 

4.1.2 Example square maps 

Consider the particular case V{x) = ^x^, and write T — Define: 

b = Vl - 12r (4-20) 



We find 



We find: 



2t 

vi^- , V2^0 , v^^-Uf (4-22) 
7 



The radius of convergence of F^^^ is |r| < ^ for all g 
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4.2 2-Matrix model 



The 2-matrix model was introduced by Kazakov [36] , as the Ising model on a random 
map. 

/ dMi dM2 e-T (Vi(Mi)+V2(M2)-Ci2 M1M2) 

where Vi(Mi) = - Vi{Mi) is a polynomial in Mi, such that Vi(Mi) contains 

only terms with degree > 3, and where V2(M2) = — V2(M2) is a polynomial in 

M2, such that V2(M2) contains only terms with degree > 3, and we assume C12 = 1: 

Indeed, it generates surfaces made of polygons of two possible colors (call them polygons 
carrying a spin + or — ) glued by their sides (no crossing here). The weight of each 
surface depends on the number of neighbouring polygons with same spin or different 
spin, which is precisely the Ising model. If Cu = C22 and Vi = V2, this is an Ising 
model without magnetic field, otherwise it is an Ising model with magnetic field. 

Let us describe the solution. 
Consider the following rational curve 



y(.)=7.-^ + E£f^'A^^ 



degy/. . (4-27) 



where all coefficients jjakyPk are determined by: 

yiz)-Vlixiz)) ~ -± + Oiz-') 

z—>oo 

x{z)-VM^)) --- + 0{z') 

(4-28) 

The coefficients 7, ak, (3k are algebraic functions of t, and we must choose the branch 
such that 7 — > at t = 0. 

The curve {x{z),y{z)) 2; G C is called the spectral curve. It is a genus zero algebraic 
curve. There are deg V2 branch points solutions of x'{z) = 0. 

Notice that we have: 

Res ydx = t = — Res ydx (4-29) 



Then one has [2]: 

= -( Res Vi{x)ydx + Res {xy — V2{y))ydx 

2 V 2— ►oo z— >0 
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-t Res Vi(x) t Res {xy - V^iv)) - In — 

2^oo ^ Z^O Z I \ ^ ) ^ 

(4 - 30) 
and [211 : 

i^^^^ = In (^i^^^ 'n ^^V^)) (4-31) 

where Oj are the solutions of x'(aj) = 0, and tdeg{V2) the leading coefficient of V2. 
The other F(s)'s are found as follows [T5] : 

Let ttj, i = 1, . . . , deg V2 be the branch points, i.e. the solutions of x'(aj) =0. If 2; 
is close to a branch-point aj, we denote ^ the other point such that 

z ai , x{z) = x{z) and ^ — > (4-32) 

Notice that z depends on the branch-point, i.e. 1 is not globally defined. We also 

define: ^ 

$(^) = / ydx (4-33) 
^ 20 

$(z) depends on the base-point zq and on the path between z and but the calcula- 
tions below don't depend on that. 

We define recursively: 

Wi^\p) = (4-34) 

Wi'\p,q) = -^ (4-35) 



(p-q) 



M^fc+\(p,Pi,...,Pfc) 

- -- y Res 7 — . . . ^ f 

2 z^ai [p — z)[p — z)[y[z) — y[z)x'[z) V 

iyif/)(z,^,Pi, . . . ,pfc) + 5^ w[%^{z,p,)w[l-^\{-z,pK/{i}) 

h=0 7c{l,2,...,fc} 

(4 - 36) 

This system is a triangular system, and all wj^^^ are well defined in a finite number of 
steps < + k. 

Then we have 1151: 



F^^') = y Res <!>{z)wi^\z) dz , g>2 (4-37) 

2 — 2q ^-^ z^ai 

The 1-matrix case is a special case of this when the curve is hyperelliptical (in that 
case z is globally defined), it corresponds to deg V2 = 2. 
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4.3 Chain of matrices 



^formal 



(4-38) 



where Vi(Mi) = - V^,(Mi) is a polynomial in Mj, such that Vi(Mj) contains only 

terms with degree > 3. The matrix C is: 



C 



-1 C22 -1 

V 



Consider the following rational curve 



1 Cpp — 1 / 



Vi = 0,...,p + 1 



(4-39) 



(4-40) 



where all coefficients jjakyPk are determined by: 



Xi+i + Xi-i = V{Xi) 

2^0 ^ 



(4-41) 



The coefficients 7, ctj^fc are algebraic functions of t, and we must choose the branch such 
that 7 ^ at t = 0. 

The curve {xi{z),X2{z)) z G C is called the spectral curve. It is a genus zero 
algebraic curve. 

Notice that we have Vz = 1, . . . ,p — 1: 



Res Xi^idxi = t = — Res Xi+idxi 

2— >00 2— >0 



(4-42) 



Then one has 



(4 - 43) 



1 / ^ 1 

-(^^ Res (Vi(xi) - -XiVi{xi))xi+idxi 

i=l 

-n. , / X 1 , ,/ / ^ ^dz , , / 7^ det C 
Res (Viixi) - -x,V:(xi)) In ^- 

i=i ^ 



F^^^ and the other F^^^'s have never been computed, but it is strongly believed that 
they should be given by the same formula as in the 2-matrix case. 
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4.4 Closed Chain of matrices 

/ dMi . . . dMp e-T ^ (^^i v>(a^.)-E?Ji MiMi+^-M^M,) 



^formal 



It is the case where the matrix C of quadratic interactions has the form: 



(4-44) 



C 



/Cn -1 -1\ 

-1 C22 -1 



(4-45) 



V-l -1 -1/ 

This model is yet unsolved, apart from very few cases {p = 2, p = 3 (Potts Model), 
p = A with cubic potentials), and there are also results in the large p limit. This model 
is still a challenge. 



4.5 0(n) Model 

JdMdM, . . .dM„e-f -(^-^-f ^ ^ 

where V contains at least cubic terms. 
We write: 

V(x) = -V{-{x + ^)) + ^{x + (4-47) 

This model can easily be analytically continued for non integer n's. Indeed, combi- 
natoricaly, it is the generating function of a loop gas living on the maps, n is the "fu- 
gacity" of loops, i.e. the n dependence of the weight of each configuration is n# loops, 
and the C dependence is (7~1^^S^^ loops Q^ie often writes: 

n = 2 cos (i/vr) (4-48) 

The 0{n) model was introduced by I. Kostov in fl3] then in [I^ [21] . and it plays a 
very important role in many areas of physics, and lately, it has started to play a special 
role in string theory, as an effective model for the check of ADS-CFT correspondence. 

The leading order 1-cut solution of this model is well known [261 123, [28] . 
For any two distinct complex numbers a and 6, define: 

m=l-- , T=—^ 4-49 
0^ K[mj 

where K{m) = K' {1 — m) are the complete elliptical integrals ([Ij). 

We also consider the following elliptical function (defined on the torus (1, r)): 

, . ., cn{2K{m)u,m) 
x(u) = ih — ; [ — f (4-50) 
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where sn and cn are the elhptical sine and cosine functions. 

Then we define the following function on the universal covering (i.e. on the complex 
plane) : 

G„M.^f(eW^^^e-W.?l^, (4.51) 

where if is a normalization constant such that: 



\imG^{u)/x{u) ^1 (4-52) 

u— >0 



It satisfies: 



G^{u+T)+G^{u-T)-nG^{u) = , G{u+1) = G{u) , G{u) = G{-t-u) 

(4-53) 

We have: 

G,{uf + G,{-uf-nG,{u)G,{-u)^{2 + n){x\u)-el) , = x(^) (4-54) 

G,_,{uf + Gi_,{-uf + nG,_,{u)G,_,{-u) = (2 - n){x^{u) - e\_,) (4-55) 

77 77 

G'i_,(«)G,(«) - G'i_,(-m)G,(-«) + -Gi_,(-«)G,(«) - -Gi_,(«)G',(-«) 
, , , 2 sin (z/vr) / , 

= -W ^n(.A-).n((l-.)A-) 

Then we define: 

^(^^) - Po.( '^"""-^-^';^-!-^--"-'"'-^-'- ) (4-57) 
.,..) . p,( P-»)(.--eU^(V^^^-MV-.C._.). ) 

where Pol means the polynomial part in x at large x (i.e. at m — >■ 0), and the subscript 
+ and — mean the even and odd part, and where 

'^^"^'^^^ (4-59) 



sn(i/i^) sn((l - i/)/s:) 



A and B are polynomials of x^ . 

Then, a and h are completely determined by the 2 conditions: 

Res (^^^^!M^^^lM± + S(a;2(ii)) G,(ii)_^ dx(ii) = (4-60) 

Res (A(x2(m)) G'i_^(m)_ + x(ii)S(x2(ii)) G'^(w)+) cix(ii) = i (4-61) 

u— »0 
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Once we have determined a, b, m, A, B, we define the resolvent: 

uj{u) = A{x{uy)Gi-^{u) + x{u)B{x{uy)G^{u) (4-62) 
which is the first term of the formal large expansion of: 

and which satisfies: 

u{u+T)+uj{u—r)+nu{u) = V'{x{u)) , u{u+l) = u{u) , uj{u) = uj{—t—u) 

(4-64) 

The free energy is then found by: 

^ (1-77) In {a^g{m)) (4-65) 



df^ ' 2 
where 

9' 1 



g m[\ — m)sn'^[uK[m)^m) 
All this is described in [30]. The special cases of n integer, and in general when u 
is rational, can be computed with more details. 

4.6 Potts model 

_ /^Mi ■■.rfMQe-T^Eti|MfH-y(M.)+gE.^,M.Af,) 
^'""^^^ " jdMdM,... dM^ e-T Et. hM^^% E?, m.m,) ^^-67) 
or equivalently: 

_ JdMdMi ...dMQe-TT>-(EaaMf+nM.)+gM^-CEa.MM.) 

The Random lattice Potts model was first introduced by V. Kazakov in [39l HQ] . 
The loop equations have been written and solved to leading order, in particular in 

[m ESI [TO]. 

4.7 3-color model 

/ dMi dM2 dM^ e-T^ C^{M^+Ml+Mi)-giM,M,Ms+M,MsA'h)) 

^'""^^ ^ / dM, dM, dM, e-T^ ihiM?+Mi+Mi)) (4-69) 

The loop equations have been written and solved to leading order, in particular in 

ESI 
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^ formal 



4.8 6- vertex model 

/ dM dM^ g-f TV (MMt-M2M2t+cos/3(MMt)2) 

^ JdMdM^e-T^(^^^^'^ 
where M is a complex matrix. The loop equations have been written and solved to 
leading order, in particular in [571 US] • 

4.9 ADE models 

Given a Dynkin diagram of A, D or E Lie algebra, and let A be its adjacency matrix 
{Aij = Aji links between i and j). We define: 

i i<j 

where Bji = B\- are complex matrices, and Mj are hermitian matrices. 

The loop equations have been written and solved to leading order, in particular in 

Halm- 

4.10 ABAB models 

/ '"l "-2 ,2 T,2 

\{dA,\{dB, ^^^^^.^^^.^^9^..^^B,A.B, (4_72) 
i=i j=i 

This model is yet unsolved, apart from few very special cases [H]. However, its solution 
would be of great interest in the understanding of Temperley Lieb algebra. 

5 Discussion 

5.1 Summary of some known results 

We list here some properties which are rather well understood. 

• The fact that formal matrix integrals are generating functions for counting dis- 
crete surfaces (also called maps) is well understood, as was explained in this 
review. 

• The fact that formal matrix integrals i.e. generating functions of maps satisfy 
Schwinger-Dyson equations is well understood. 

• The fact that formal matrix integrals and convergent integrals don't coincide in 
general is well understood. In the examples of the 1-matrix model, 2-matrix Ising 
model or chain of matrices, it is understood that they may coincide only in the 
"1-cut case" , i.e. if the classical spectral curve has genus zero. 
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The fact that the formal 1-matrix, 2-matrix or chain of matrices integrals are r 
functions of some integrable hierarchies is well understood too. 

For the formal 1-matrix model and 2-matrix Ising model, all F*^^)'s have been 
computed explicitly. The result is written in terms of residues of rational func- 
tions. 

For the chain of matrices, F^^^ is known explicitly |25j, and it is strongly believed 
that all F^9)^s are given by the same expression as for the 2-matrix model. 

Multi-cut formal matrix models are well studied too, and they can be rewritten 
in terms of multi-matrix models. For the 1 and 2 matrix models, the expressions 
of the F^3)^s are known explicitly (in terms of residues of meromorphic forms on 
higher genus spectral curves). 

Some open problems 

the large N asymptotics of convergent matrix integrals, are not directly relevant 
for combinatorics, but are a very important question for other applications, in 
particular bi-orthogonal polynomials, universal statistical properties of eigenval- 
ues of large random matrices, and many other applications. This question is 
highly non-trivial, and so far, it has been solved only for the 1-matrix model 
[6], [19], [22], and a few other cases (e.g. the 6- vertex model [8]). The method 
mostly used to prove the asymptotic formulae is the Riemann-Hilbert method 
[HI [I9j, which consists in finding a Riemann-Hilbert problem for the actual ma- 
trix integral and for its conjectured asymptotics, and compare both. There is a 
hope that probabilists' methods like large deviations could be at leats as efficient. 

For the 2-matrix model and chain of matrices, the topological expansion of formal 
"mixed" expectation values (e.v. of invariant monomials whose words contain at 
least 2 different letters) has not yet been computed. This problem is a challenge 
in itself, and has applications to the understanding of "boundary conformal field 
theory". In terms of combinatorics, it corresponds to find generating functions 
for open surfaces with boundaries of prescribed colors. 

Many matrix models have not yet been solved, even to planar order. For instance 
the closed chain of matrices where Cij = + + CtSa and Cpi = Cip = 1. 
For the Potts model, 6- vertex model, 3-color model, 0(n) model, ADE mod- 
els, only the planar resolvent is known. For the AnBmAnBm, almost nothing is 
known, although this model describes the combinatorics of Temperly Lieb alge- 
bra. 
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• Limits of critical points are still to be understood and classified. Only a small 
number of critical points have been studied so far. They have been related to 
KdV or KP hierarchies. Critical points, i.e. radius of convergence of the t-formal 
power series, are in relationship with asymptotics numbers of large maps (through 
Borel's transform), and thus critical points allow to count maps made of a very 
large number of polygons, i.e. they can be interpreted as counting continuous 
surfaces (this was the idea of 2D quantum gravity in the 80's and early 90's). 
This is yet to be better understood [16], U\- 

• extensions to other types of matrices (non-hermitian) have been very little studied 
compared to the hermitian case, and much is still to be understood. For instance 
real symmetric matrices or quaternion-self-dual matrices have been studied from 
the begining |38], and they count non-orientable maps. Complex matrices, and 
normal complex matrices have played an increasing role in physics, because of 
their relationship with Laplacian growth problems [54], or some limits of string 
theory [1]. Complex matrices count maps with arrows on the edges (see the 6- 
vertex model [SZ]). Other types of matrix ensembles have been introduced in 
relationship with symmetric spaces [58], and it is not clear what they count. 

• And there are so many applications of random matrices to physics, mathematics, 
biology, economics, to be investigated... 
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